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^ Abstract 

fSJ .We study a general problem how the gap in a nonmagnetic band insulator closes by tuning a parameter We review our recent results 
^jQbn the classification of all the possible gap closing in two and three dimensions. We show that they accompany the change of Z2 
^ topological numbers, and that the gap closings correspond to phase transitions between the quantum spin Hall and the insulator 
phases. Interestingly, in inversion-asymmetric three-dimensional systems there appears a gapless phase between the quantum spin 
Hall and insulator phases. This gapless phase is due to a topological nature of gap-closing points in three dimensions, but not in 



two dimensions. 
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1. Introduction 

The physics of pure spin current has been attracting much at- 
tention recently in various fields of condensed matter. Pure spin 
fcurrent is a flow of electron spins which is not accompanied by 
the charge current. As is difi'erent from the spin itself, spin cur- 
rent is even under the time-reversal operation, and thus can be 
nonzero in nonmagnetic systems. Thereby the spin current has 
bpened up new fields in nonmagnetic condensed materials. One 
example is the intrinsic spin Hall effect (SHE) \X, 2] driven by 
the spin-orbit interaction. The SHE has offered us a new way to 
inject and detect spin currents in metals and semiconductors. 

The physics of topological insulators, also called the quan- 
tum spin Hall (QSH) systems , is another example of 
new physics due to the pure spin current. Two-dimensional 
(2D) QSH systems are insulators in the bulk, while they have 
gapless edge states. Similarly, three-dimensional (3D) QSH 
systems (topological insulators) are insulators with gapless sur- 
face states. One of the important and novel aspects of these sys- 
tems is that these gapless edge/surface states are robust against 
perturbation which respects time-reversal symmetry. Namely, 
they remaingapless even when nonmagnetic impurities or dis- 
order exist iaBlj which causes novel transport properties JMSl- 
The robustness of these gapless states can be interpreted as pro- 
tected by topology. The QSH systems are characterized by non- 
trivial Z2 topological number The Z2 topological number takes 
only two dififerent values for v: v = (also called v =even) 
or V = 1 (also called v =odd). When v = 1, the system is in 
the QSH phase, while when v = the system is in the insulat- 
ing phase. This topological number remains unchanged as long 
as the bulk gap is open and the time-reversal symmetry is pre- 
served. Transitions between ordinary insulator (v - 0) and the 
QSH system (v = 1) occur only when the gap is closed. The 
topological classification has been discussed in various theo- 
retical papers, but they might not be easily accessible for non- 
experts in topology. 



Instead of directly dealing with the Z2 topological num- 
bers, we take a different route and study the physics of gap- 
closing. This has been studied in a series of papers of the au- 
thor 1 U, 121. In the present paper, we summarize these 



results in these papers. In generic nonmagnetic band insulators, 
we examine whether the gap closes when we vary a parame- 
ter in the system. As a result we can classify all the possible 
types of gap-closing. In fact, in every case of the gap closing, 
we show in this paper that the Z2 topological number changes. 
Namely, the gap-closing physics is equivalent to the phase tran- 
sition between the QSH system and the ordinary insulator. As 
a byproduct, we obtain effective theories focusing on the vicin- 
ity of the gap-closing points at the transition. This gap-closing 
physics involves local features in k space, and is equivalent to 
global topological structure in k space , 13 ,[illl5lil6tl. We 
note that this gap-closing physics tells us only the change of 
the topological number, and we cannot tell which side of the 
gap-closing is the topological insulator and which side is thee 
trivial insulator. This is the characteristic of the theories involv- 
ing local features in k space. On the other hand, the gap-closing 
physics gives us a number of new results, such as effective the- 
ories expanded in terms of k, or the universal phase diagrams, 
as we discuss in this paper 

As a result of our analysis, we get a full understandingof the 
gap-closing physics, which is different between 2D flUS] and 
3D [Ts*, To], and between inversion-symmetric and inversion- 
asymmetric systems. The most interesting is the 3D inversion- 
asymmetric systems; we find that there necessarily exists a fi- 
nite region of gapless phase between the QSH phase and the 
ordinary insulator phase. This gapless phase is of topological 
origin, as we describe in this paper. Henceforth we consider 
only clean systems without any impurities or disorder The 
time-reversal symmetry is assumed throughout the paper. In the 
present paper we only deal with time-reversal-invariant systems 
with spin-orbit coupling, without assuming additional symme- 
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tries such as spin-rotational symmetry. Therefore, the systems 
without spin-orbit cou plin g, such as the organic system with 
bulk Dirac-like bands [17], belong to a different class of sys- 
tems (orthogonal class) and are beyond the scope of the present 
paper. 

2. Gap closing at the phase transition 

Putting topological insulators aside for a while, we ask our- 
selves a general question when and how the gap of a nonmag- 
netic band insulator (with spin-orbit coupling) closes by tuning 
a single parameter m in the system. We assume that the system 
is time-reversal symmetric. Depending on the system consid- 
ered, this parameter m can have any physical meaning, such as 
pressure, chemical content, or interatomic distance. Because 
our goal is to pursue the phase transition, we will focus on gap- 
closing achieved by tuning a single parameter m in generic sys- 
tems. We call such gap closing as "generic gap closing". The 
word "generic" means that any terms allowed by symmetry are 
set to be nonzero. There are several reasons why we focus on 
generic gap closing. First, in real materials there are many fac- 
tors which determine the band structure. Therefore, any kind 
of terms in the Hamiltonian are set to be nonzero, as long as 
the terms are allowed by symmetry. In other words, we assume 
that there are no accidental degeneracies unless it is required 
by symmetry. Second, only the generic gap closings are nec- 
essarily related with phase transitions. Nongeneric ones do not 
correspond to phase transition, because such kind of nongeneric 
gap-closing may disappear by small perturbation. (If such non- 
generic gap closing does not disappear by perturbation, it is 
nothing but a generic gap closing studied in this paper. Thus 
it is enough to study only the generic gap closing.) By these 
reasons, non-generic gap-closing which requires fine tuning of 
some other parameters is excluded in our analysis. 

With these reasonable and general assumptions, we consider 
a problem whether the gap will close by tuning only a single 
parameter m in the system. In respective cases considered in 
the following, we find necessary conditions for the gap to close, 
and see whether it can be satisfied by changing a parameter m. 
We consider a gapped spin- 1/2 time-reversal-symmetric system 
with spin-orbit interaction. A Hamiltonian matrix for Bloch 
wavefunctions can be written in a block form. 



i/(k) = 



(1) 



/2iT(k) /2u(k) 

The dimensions of the Hamiltonian i/(k) can be arbitrary for 
the subsequent analysis, but for simplicity we set the dimen- 
sion to be minimal, as long as the gap-closing physics is ap- 
propriately captured; we retain only the states which are in- 
volved in the gap closing. In the following, we find that the 
dimension turns out to be four (see ^2.11 ) and two (see ^2.21 ) 
for the inversion-symmetric and inversion-asymmetric systems, 
respectively. The spectrum is assumed to have a gap, within 
which the Fermi energy lies. Because the time-reversal opera- 
tor is given by = io-yK with K being complex conjugation, 
the time-reversal-symmetry results in 



which is cast into the equations /!-r|(k) - /z^(-k) and /i-rx(k) = 

Before proceeding to the specific cases, we introduce the 
momenta with k = -k (mod G) where G is a reciprocal lat- 
tice vector Such momenta are called the time-reversal invari- 
ant momenta (TRIM) F, 

lie HE a 

and have the val- 
ues r,=(„,„j„,) = (nibi + n2b2 + n3b3)/2 in 3D, and r,=(,„„,) = 
(«ibi + «2b2)/2 in 2D, where nj = 0, 1 and bj are primitive 
reciprocal lattice vectors. These momenta are invariant under 
time-reversal, and are important in the subsequent discussion. 

Before proceeding to the cases with inversion-symmetric and 
inversion-asymmetric systems, we explain the definition of in- 
version symmetry. The space inversion means that the sign of 
every space coordinate is changed, and the crucial point here is 
that the inversion leaves spins unchanged. In three dimensions 
it means x — > -x, y — > -y, and z — > -Z- In two dimensions 
the space coordinates change as jc — > -x, y — » -y, while it is 
different from the rotation within xy plane, because the space 
rotation changes — » -s-\ s-' — > -s', and s~ — > s'", while in- 
version does not change spin. To see this more clearly, we take 
the Rashba spin-orbit term A{s x k), as an example, where s is 
the spin. This term changes sign under inversion since k flips 
its sign while s does not. This sign change under inversion can 
be regarded as a reversal of the +z direction, which represents 
structural inversion asymmetry causing the Rashba coupling. 

2.1. Inversion symmetric systems 

In the inversion-symmetric systems, the energies are dou- 
bly degenerate at each k. The gap closing involves a doubly- 
degenerate valence band and a doubly-degenerate conduction 
band. Hence, four states are involved in gap closing, and we 
consider 4x4 Hamiltonian matrix //(k) (/i(j,^(k) in Eq. ([T]) is 
2x2). The inversion-symmetry is expressed as 



H(-k) = PH(k)p-\ M(-k) = Pu(k), 



(3) 



where Pisa unitary matrix independent of k, and M(k) is the pe- 
riodic part of the Bloch wavefunction: ^k(r) = u{k)e'^''^. As the 
inversion does not change spin, this unitary matrix P is block- 
diagonal in spin space: 



P = 



(4) 



By using O) and after judicious unitary transformation, all 
cases are shown to reduce to a case Pi = Pi - diag(77fl, rjh) with 
T]a = +1, rji, - +1. rja and rjh represent the parity eigenvalues of 
the atomic orbitals involved. 

Occurrence of gap closing turns out to be different for the 
cases (i) T]a - rjb and (ii) t/a = -rib. The case (i) rja - rjt - 
+ 1 means that the atomic orbitals a, b have the same parity, 
such as two 5-like orbitals or two /?-like orbitals. The inversion 
symmetry (O then imposes that 



/i^ is an even function 



of k, and h^^ - is an antisymmetric matrix, even function 
of k. The generic Hamiltonian becomes 



i/(k) = o-yW (-k)cr,, 



(2) 



i/(k) = £o(k) + ^fl,(k)r,- 



a-i a\ — ia2 —04 - ia^ ^ 

ai + ia2 —03 04 + ifls 

04 - ias fli + /a2 

— 04 + /as fli— ;fl2 -fl3 



where fl,'s and £0 are real even functions of k. F, are 4x4 matri- 
ces forming Clifford algebra, given by Fi = 1 ® Tv, F2 - o-,®Ty, 
Fj = 1 18> Tr, F4 = cTy ig) Ty, aiid F5 = cr J. ig) T,., where cr; and t/ are 
Pauli matrices acting on spin and orbital spaces, respectively. 

The eigenenergies are given by £0 ± ^Yj^=i <^j- The g^P closes 
when fl,(k) - for i - 1, ■ ■ ■ , 5. Namely, these five independent 
equations should be satisfied for the gap to close. This number 
is called codimension. We see that in general there are no solu- 
tions of kjt, ky and m satisfying these five conditions. It follows 
from our assumption of "generic" systems, i.e. we exclude a 
possibility that these five conditions are satisfied by accident by 
only three available parameters k^, ky and m. Hence, the gap 
never closes for (i) 77„ =77^. This is caused by the level repul- 
sion between the states with the same parity. As we can see 
from this example, when the codimension exceeds the number 
of available parameters, we cannot expect the gap closing to oc- 
cur in generic systems. Thus the codimension represents how 
hard it is to achieve gap closing, by overcoming the level repul- 
sion between the conduction and the valence bands. 

The result is diff'erent for the other case (ii) rja = -rji, - +1, 
i.e. P - TjgT, - +T,, where the two atomic orbitals have differ- 
ent parity. After a simple calculation, one can write the Hamil- 
tonian as 



i/(k) = flo(k) + a5(k)F5 + 



(5) 



where ao(k) and a5(k) are even functions of k, and /7*'*(k) are 
odd functions of k. Here FJ are 4x4 matrices forming Cliff'ord 
algebra, given by F'j = cr. ® t^^T'^ - 1 ® t,,, F^ - cr^ (gi t^, 
F4 = CTy ® Tf, and F^ = 1 ® t~. In this case the gap closes when 
five equations fl5(k) = 0, b^'\k) = are satisfied. At generic 
k, these five equations cannot be satisfied simultaneously even 
when a single parameter m is changed, again by our assumption 
of generic systems. On the other hand, at the TRIMs k = F,-, the 
situation is different. At these points k = -k (modG) holds, and 
the odd functions b^'\k) vanish identically, and one has only to 
tune a5(k) to be zero. Thus, the gap closes at k = F, by tuning 
a single parameter m. By putting Ak = k - F,-, the Hamiltonian 
is expanded to linear order in Ak as 



H(k) ~Eo + mr',+Y, [l^'^ ■ Ak) F;, 



(6) 



/■=! 



where £0 and m are constants, and yS*'* (/ = 1 , ■ ■ ■ , 4) are two- 
dimensional real constant vectors. In this model, m is the pa- 
rameter which controls the gap closing, and the gap closes at the 
TRIM F,- when m = 0. After judicious unitary transformations, 
the Hamiltonian finally becomes block-diagonal; 



H{k) = £0 



m 



Z- 

—m 



m 

-Z- 



-z+ 
-m 



(7) 



where z± - biAk^ + bj,Aky + ib2^ky with real constants bi, 
b2 and b^. Note that in materials with e.g. 3- or 4-fold 
rotational symmetry within xy plane, one has bi - b2 and 
^3 - 0, leading to z± Ak^ + iAky. Thus we have shown that 
the generic Hamiltonian with spin-orbit coupling with time- 
reversal- and inversion-symmetries decouples into a pair of 
Hamiltonians describing two-component fermions with oppo- 
site sign of the corresponding mass terms. The eigenenergies 
are E = Eq ± yjm?- + z+Z-. The inversion matrix in this basis is 
written as rja - ?7ndiag(l, -1, 1,-1). 

In 3D, by the similar analysis, we can write down the generic 
form of the effective Hamiltonian. If we assume an additional 
symmetry of e.g. 3- or 4-fold rotational symmetry within xy 
plane, the Hamiltonian simplifies to 



H(k) = £0 + 



m 



CAL \ 



Z- 

—m CAkr 

CAk, m — z+ 

, CAk. -z- -m 



(8) 



where z± - A{Ak^ ± iAky) with real constants A and C, after 
unitary transformations. By a unitary transformation, this is the 
same as the Hamiltonian for Bi2Se3 discussed in fis'] up to the 
linear order in Ak. This model is used in |19, 2Q] to discuss 
disorder effects on the 3D topological insulators. 

2.2. Inversion asymmetric systems 

In the inversion-asymmetric systems, the phase transition is 
quite different from the inversion-symmetric systems. First we 
show that the gap cannot close at the TRIMs k - F,. This 
is because at the TRIMs, all the states are doubly degenerate, 
and the codimension is five |21 , 22 ] . To see this explicitly, we 
consider 4x4 Hamiltonian matrix with the constraint (|2]i, as the 
number of states involved is four. It leads to a result 



i/(k = F,) = £o + ^fl,r; 



(9) 



where a,'s and Eq are real. Its eigenenergies are given by Eq ± 

-^2^=1 '^r The gap between the two (doubly-degenerate) bands 
close when a, = for / - 1, ■ ■ ■ , 5, but they cannot be satisfied 
simultaneously by tuning only one parameter m. (We note that 
k is fixed to be k = F,.) Therefore, the gap does not close at the 
TRIMs k = F,. 

If k is not among the TRIMs, each band is non-degenerate. 
We thus retain only one valence and one conduction bands 
(Therefore the states are involved and i/(k) is set to be 2x2). 
Equation ^ gives no constraint for the Hamiltonian near k. 
The Hamiltonian near k thus reduces to 



H = 



a 

c* 



(10) 



where a, b are real functions of k and m, and c is a complex 
function of k and m. The gap closes when the three conditions 



a — b, Rec = 0, Imc = 0, 
are satisfied, i.e. the codimension is 3 j23, 24]. 



(11) 
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Figure 1: Generic gap-closing in 2D for (a) inversion-asymmetric and (b) 
inversion-symmetric cases. In case (b) all the states are doubly degenerate by 
Kramers theorem. 



In 2D there are three variables k^-, ky and m, which is the 
same number as the gap-closing conditions; therefore, the gap 
can close at some F,) when the parameter m is tuned. Let 
mo to be the value of m at the gap closing. Because of the 
time-reversal-symmetry, the gap closes simultaneously at k = 
±ko, at m - niQ, as depicted in Fig. [T](a). In 2D, the effective 
Hamiltonian near the gap-closing point k = ko(9i F,) reduces to 

'H = Eo(m,k^-,ky)+(m-mo)o-,+(k^-kM)o-^+(ky-kyo)o-y(12) 

after unitary and scale transformations ifioll . Thus, to summa- 
rize our theory in the 2D systems, generic gap-closings are 
classified into two cases shown in Fig. [T] (a)(b) for inversion- 
asymmetric and inversion-symmetric cases respectively. Other 
types of gap-closings are prohibited because of level repulsion 
between the states in the valence and the conduction bands by 
tuning a single parameter m. 

In 3D, in contrast with the 2D case, there are 4 variables 
m, kjc, ky, k~, and the situation is different. Because the gap- 
closing condition consists of three equations, it determines a 
curve in the four-dimensional (k^, ky, k^, m)-hyperspace, which 
we call a "string" C. Generally, this string C in (k^, ky, k,, m)- 
hyperspace occupies a finite region in m-direction mi < m < 
Within this region m\ < m < 1112 the bulk gap is closed and 
the system is in a gapless phase. As we see in the subsequent 
discussion, this gapless phase is a topological phase, protected 
by the existence of monopoles and antimonopoles in k space. 

2.2.1. Monopole-antimonopole pair creation and annihilation 
in k space 

We note that we are focusing on the 3D inversion-asymmetric 
systems. To see the behavior of the gap-closing points in k 
space, we note that each gap-closing point carries a topologi- 
cal number. Such gap-closing point is regarded as a monopole 
or an antimonopole in k space 1 25 , 26 , 27 , 28 , 2^ . It can be 
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m=m„ 



k=r: 



Figure 2: Trajectory of the gap-closing points for (a) inversion-asymmetric and 
(b) symmetric systems. For (b) inversion-symmetric systems, the gap-closing 
point is located at k = F,-, and is an isolated point in the m-k space. Only 
at m = mo the system is gapless. For (a) inversion-asymmetiic systems, on 
the other hand, the gapless points are created in monopole-antimonopole pairs 
at m = mi, and move in k-space as m is varied. Solid and broken curves 
represent the trajectories for monopoles and antimonopoles, respectively. These 
two curves together form a closed loop, which is the string C in ^2.21 The 
system opens a gap only by pair annihilation of these gapless points at m = m2. 



seen by associating the Bloch wavefunctions with nondegener- 
ate spectrum to U(l) gauge field in the k space; 



A„(k) = 
B„(k) = 

P«(k) = 



-'XMk«|Vk|Mk«), 

Vk X A„(k), 



1 

2^ 



Vk ■ B„(k). 



(13) 
(14) 

(15) 



where |Mk,n) is the Bloch wavefunction for an n-th band. A„(k), 
B„(k) are called Berry connection and Berry curvature, respec- 
tively. p„(k) is called a monopole density. p„(k) vanishes 
when the «-th band is not degenerate with other bands, because 
Vk ■ (VkX ) = 0. On the other hand, p„(k) does not vanish when 
the n-th band touches with another band in energy at some k- 
point k - k(). In such case, the Bloch wavefunction cannot be 
expressed as a single continuous function around k - ko; k 
space should be patched with more than one continuous wave- 
functions Isoll . as is similar to the vector potential around the 
Dirac monopole 13111 . As a result, Pn(k) is shown to have a 
^-function singularity at at k = ko S S Ezt IM llSt] . For 
example, for the gap closing at k = ko with linear dispersion 
(Weyl fermion) 



'H = £o(k) + ^/i(k)cr,- 



(16) 



where /j(k - ko) - 0, the monopole density for the lower 
band is given by ^(5(k - ko), where q = sgn(det(||),j)|k=ko 
(= +1) is called monopole charge. In general the monopole 
density has the form p(k) = 2/ qi6{\i - k/) where the monopole 
charge qi is quantized to be an integer. It is shown that when 
we vary the system by changing a parameter continuously, the 
monopole charge is conserved, because of the quantization of 
the monopole charge. The only chance for the monopole charge 
to change is a pair creation or a pair annihilation of a monopole 
(qi - I) and an antimonopole (qr = -1). 

In the present problem of nonmagnetic insulators, because 
of the time-reversal-symmetry, the distribution of monopole 



4 



charges is symmetric with respect to k = F,-: pQ.(k) - pa.(-k) - 
PQ,(2r, - k), where a is the label which is a time-reversed la- 
bel from a. From these arguments we can see that the sim- 
plest case of the gap closing is as shown in Fig.|2fa) lllllll2ll . 
Two monopole-antimonopole pairs are created at k = +ko + F, 
(ko 0) simultaneously when m = mi, and the system be- 
comes gapless. When m is increased further, the monopoles 
and the antimonopoles move in the k space, while the distri- 
bution of the monopole charges remains symmetric with re- 
spect to r,. This system can open a gap again only when 
all the monopole and antimonopole annihilate in pairs. This 
occurs at m = m2 as shown in Fig. [2] (a). Thus the over- 
all feature of the phase transition is schematically expressed 
as in Fig. [3] As the monopoles and antimonopoles are gap- 
closing points in k space, the trajectory of the monopoles and 
antimonopoles is nothing but the string C (see ^2.21 ). describ- 
ing the set of parameters (k^, ky, k,, m) which satisfy the gap- 
closing conditions ( fTTT ). Because the monopole charge is con- 
served, the monopoles and antimonopoles are created and an- 
nihilated only in pairs, which means that the trajectory C of 
the monopoles and antimonopoles form a closed loop in the 
(k, m) space. Namely, the string C has no end point, because 
an end point of C would violate the conservation of monopole 
charge. When nii < m < mi the system is gapless because 
there are monopoles and antimonopoles, which are gap-closing 
points (see Fig.|3a)). 



(a) 



/Iv 







m<mi 
(b) 



1 




1 


k=r/ * 

\ 




m<mo 



m=mo 



m>mo m 



Figure 3: Phase transition in 3D between the quantum spin Hall (QSH) and in- 
sulating phases for (a) inversion-asymmetric and (b) inversion-symmetric cases. 
In the case (b) all the states are doubly degenerate. 



3. Change of the Z2 topological number at the gap-closing 
point 

In §2 we classify all the generic gap-closings by tuning a 
single parameter in nonmagnetic insulators. In this section we 
relate this to topological insulators. We will see that all these 
gap-closings accompany the change of the Z2 topological num- 
ber and thus describe phase transition between the QSH phase 
and the ordinary insulator phase. As a result the universal phase 



diagram is as shown in Fig. |4] which is the main result of the 
paper 

3.1. Z2 topological numbers: review 

We first review the exp ressions for the Z2 topological num- 
ber given in il3l Il4 11511 . both in 2D and in 3D. We assume 
that the spectrum of the Hamiltonian has a gap, within which 
the Fermi energy Ef is located. From the ICramers theorem 
for time-reversal-symmetric systems, ;/ inversion-symmetry is 
broken, double degeneracy occurs only at the four TRIMs, and 
non-degenerate at other points. In such systems, the Z2 topo- 
logical number v is determined as 



Vdet[w(r,)] 

pf[w(r,)] 



±1. 



(17) 



Here w(k) is a unitary matrix with elements given by w„,„(k) = 
(M_k,m|0|Mk,«), and |Mk,«) is the Bloch wavefunction of an ra- 
th band whose eigenenergy lies below Ef. The branch of the 
square root of the determinant is so chosen that the wavefunc- 
tions (including their phases) are continuous in the whole Bril- 
louin zone. On the other hand, in inversion-symmetric systems, 
the formula simplifies drastically; it is given by 



(-D^^fj^,, 5, = f]^2,„(r,), 



(18) 



m=l 



where ^imiTi) (- +1) is the parity eigenvalue of the Kramers 
pairs at the TRIM F,, and is the number of Kramers pairs 
below Ef. 

In 3D, there are four Z2 topological numbers written as 
vo; (V1V2V3) OlSllla], given by 



(-1)"° -f]^,-, (-1)" 



|~~[ ^i=(ni«2"3)- 



(19) 



When V() = it is called the weak topological insulator (WTI), 
while when vq = 1 it is called the strong topological insulator 
(STI). These topological numbers in 3D determine the topology 
of the surface states for arbitrary crystal directions |15]. We 
note that among the four Z2 topological numbers in 3D, only vq 
is robust against nonmagnetic disorder, while vi, V2, and V3 are 
meaningful only for a relatively clean sample ifisll . 

3.2. Change of the Z2 topological number at gap closing 

We have shown in lflol,[Tll[l2ll that all the types of gap closing 
found in §2 accompany the change of the Z2 topological num- 
ber, and thus entail phase transition between the QSH phase and 
the ordinary insulator phase. Following 1 10. LI. lZl we explain 
its proof and its implications. 

For the inversion-symmetric systems in 2D and 3D, from 
Eq. ( fT9] l. the Z2 topological numbers are given by the parity 
eigenvalues of the occupied states. The gap at k = F,- collapses 
when m - 0. Hence only the parity eigenvalue at k = F, can 
change at the phase transition. Since the inversion matrix is 
given by P = 77fl 181 t~ - rjacro ® r,, the parity eigenvalues of (jT) 
at k = F,- are -riai- rjt,) and +770 for the lower-band states at 
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(a) 3D 



(b)2D 



QSH 




Figure 4: Universal pliase diagram for the quantum spin Hall (QSH) and in- 
sulator (I) phases in (a) 3D and (b) 2D. m is a parameter driving the phase 
transition, and S represents inversion symmetry breaking. 5 = corresponds to 
the inversion-symmetric system. 



m > and m < 0, respectively. Hence, the parity eigenvalue 
changes sign, and the Z2 topological number v changes by one. 
Thus, on the two sides of the gap closing, m > and m < 0, 
one of the phases is the QSH phase, while the other one is the 
ordinary insulator phase. 

For the inversion-asymmetric 2D systems, the homotopy 
characterization in Ref. fl6^ is applicable; for the lower band 
at the critical value m = otq, there is one vortex at ko and one 
antivortex at -ko. Thus, when the parameter m is tuned across 
m = niQ, the Chern number for the whole contracted surface 
lfl6ll changes by one. Thus, the Z2 topological numbers are dif- 
ferent by one for the m > mo and the m < mo sides. One of 
the phases is the QSH phase, while the other is the ordinary 
insulator 

For the inversion-asymmetric 3D case, we can relate the 
shape of the trajectory ("loop") of the gap closing points 
(Fig. |2la)) with the change in the topological number It has 
been discussed in II12II . The rule is simple; if the loop winds 
around a TRIM F,- once (as in Fig. |2la)), the index <5, changes 
between the phase in m > m2 and that in m < m\. The change 
of the Z2 topological numbers vq; (vi V2V3) follows from it. In 
particular, vq changes when the loop winds around only one of 
the TRIMs as in Fig.|2ja). It is reasonable since by restoring the 
inversion symmetry via continuous change of the Hamiltonian, 
Fig. 12 a) should reduce to Fig. |2b). 

To summarize the universal phase diagram is obtained as 
Fig.|4] We note that for inversion-asymmetric 3D systems there 
should be a finite region of "topological" gapless phase between 
the topological and trivial phases. This is exactly the region 
where the string C lies (see Fig.|2la)). 



4. Examples of the gap closing: models and materials 

4.1. 2D systems 

The Kane-Mele model on the honeycomb lattice is stud- 
ied within our theory and it agrees with our results. The hop- 
ping t term and the spin-orbit A^o term preserves the inversion 
symmetry. Nevertheless, the phase diagram (Fig. 1 in fs*]) is 
in the A^-Ar space, and both /I,, and terms break inversion 
symmetry. Thus this model in Ay-A^ space falls into the class 
of inversion-asymmetric 2D model (Fig. [T] (a)). As we have 
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-0.2 -0.1 0.1 0.2 -0.1 -0.05 0.05 0.1 

Figure 5: Phase diagrams for the Fu-Kane-Mele model with St^ = 0, St^ = 0. 
t\ and t2 are the bonds along the [111] and [111] directions. We put /iso = 0. If. 
The axes are in the unit of t. (a) The phase diagram in St\ -St2 plane obtained in 
Ref. [15]. Ay is set as zero, (b) The phase diagram in the 6t+-Ay plane. We have 
introduced A,, into the Fu-Kane-Mele model. Here St+ = dt[ + St2 is changed, 
while we fix St- = Sti - Sti =0.1?. The arrows in (a) and (b) correspond to the 
identical change in parameters. 



predicted, the gap then should close at non-TRI]VI points. In ac- 
cordance with this expectation, the gap closes at the K and K' 
points (i.e. non-TRIM) at the phase transition. 

In the CdTe/HgTe/CdTe quantum well, the 2D quantum spin 
Hall state has been observed experimentally I32i, i8il. By increas- 
ing the well thickness d over the critical thickness dc ~ 64A the 
system undergoes a transition from the ordinary insulator to the 
QSH system. The gap closing here is the class of 2D inversion 
symmetric system in Fig. [Hb), and should occur at one of the 
TRIMs. Indeed, when d - dc, the bulk gap closes at the F point 
(i.e. TRIM), and the parity eigenvalues are exchanged. More- 
over, the effective model is discussed in ll32ll up to the order. 
We can compare it to our effective theory (|7]) and they are in 
complete agreement to the linear order in k. 

4.2. 3D systems 

We take the 3D tight-binding model proposed by Fu, Kane 
and Mele Iil5i] on a diamond lattice as an example. This model 
shows a transition between STI and WTI. This model is 3D 
inversion symmetric, and the gap should close at one of the 
TRIM. Indeed, as calculated in [13,1 , the band gap closes at one 
of the X points. Furthermore, one can use this model to confirm 
the universal phase diagram in Fig. ID including the topological 
gapless phase for inversion-asymmetric 3D systems. We ex- 
plain it, following the analysis in Ref |[l2ll . We break inversion 
symmetry by adding a staggered on-site potential Ay, and cal- 
culate the phase diagram. The results are shown in Fig. |3b), 
where the vertical axis represents the staggered on-site poten- 
tial, corresponding to the parameter S in Fig.lUa), Thus we can 
see that if we break the inversion symmetry there occurs a fi- 
nite region of a gapless phase in the phase diagram. These are 
in accordance with our theory. Furthermore, The trajectory of 
the gap closing points (not shown in this paper, details are in 
lfl2ll ') agrees with Fig. |2ja). Thus we have seen that the uni- 
versal phase diagram in Fig.|4]holds true in the Fu-Kane-Mele 
model. 

Bi2Se3 and Bi2Te3 El I33I1 are topological insulators un- 
der intensive experimental and theoretical research. Their bulk 
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bands have a direct gap at the F point. The highest valence 
band and the lowest conduction band have the opposite parities 
at the r point, and if the spin-orbit coupling is weakened their 
energies will be reversed and they become ordinary insulators. 
In other words, topological insulator phases in these materials 
are due to the gap closing at the F point by making the spin- 
orbit coupling larger These are inversion-symmetric, and are 
classified to the 3D inversion-symmetric systems, whose gap 
closing occurs only at the TRIMs, in agreement with the above 
discussion. The eff'ective model is constructed in Iil8.] . and is in 
agreement with our theory in ^ u p to linear order in k. 

Both in 2D and in 3D, bismuth 134,13511 . antimony and their 
alloys 1361 also are interesting materials for topological insu- 
lators luHo. The 3D Bi has v =0, while the 3D Sb has v =1. 
This is opposite to the case of the 2D fsTj. (Note that although 
Bi and Sb are semimetals, there are direct gap everywhere in 
the Brillouin zone. This makes Z2 topological invariants well- 
defined, by artificially making the bands to be fully gapped by 
continuous deformation of the bands to eliminate the small car- 
rier pockets.) To resolve this, numerical analysis was performed 
in [37] by artificially changing the interlayer hopping weakened 
by a factor / (0 < / < 1). Between / - 1 (3D limit) and / = 
(2D limit), the gap closes at several parameter values of /, and 
this is described by the gap closing in inversion-symmetric 3D 
systems. The gap closes only at the TRIMs, in accordance with 
our theory. 

5. Conclusions and Discussions 

To summarize, we analyze gap closings in 2D and 3D non- 
magnetic insulators by tuning a parameter We obtain the uni- 
versal phase diagrams in 2D and 3D, as schematically shown 
in Fig.m in a plane of the control parameter m and another pa- 
rameter S representing an inversion-symmetry breaking. The 
overall behavior of gap closing is schematically shown in Fig.[T] 
for2DandFig.|3]for3D. 

We find that in inversion-asymmetric 3D systems, there lies 
a finite region of the gapless phase between the quantum spin 
Hall and ordinary insulator phases. We checked this for the 
Fu-Kane-Mele model. We described the phase transition in 
terms of the motion of the gap-closing points (i.e. monopoles 
and antimonopoles) in k space. The gapless phase in the 
inversion-asymmetric 3D system originates from the conserva- 
tion of "monopole charge". These cases of gap-closing exactly 
coincide with the phase transitions between the QSH and the in- 
sulating phases. In this sense our theory characterizes the QSH 
phase from the local features in k space. All the known models 
exhibiting phase transition between the two phases are special 
cases of this general classification. To confirm this scenario, 
experimental pursuit of 3D inversion-asymmetric topological 
insulators will be interesting and important. 
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